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1. Introduction 

Matrix theory [1J is a remarkable proposal for a non-perturbative definition of M- 
theory, a supposedly consistent 11-dimensional quantum theory having supergravity as its 
low energy limit. Matrix theory is defined as the maximally supersymmetric quantum 
mechanics of U(N) matrices, describing the lowest states of open strings connecting N 
DO-branes. Supergravitons, for instance, appear as bound states of these DO-branes. The 
original conjecture of Banks, Fischler, Shenker and Susskind pj relates the large iV limit of 
Matrix theory to M-theory in the Infinite Momemtum Frame. Later, Susskind expanded 
the conjecture to a relation between finite iV Matrix theory and M-theory with a compact 
null-direction. 

There have been several successful tests of these conjectures so far, see, e.g., for an 
overview. Part of the original evidence involved comparing graviton scattering phase shifts 
computed in Matrix theory and supergravity in the limit of large impact parameter. The 
fact that these agree originates from the work of Douglas, Kabat, Pouliot, and Shenker 
0, where it was shown how supersymmetry leads to a cancellation of the contribution 
from the massive modes of the strings connecting the DO-branes, implying that the long 
distance behavior can be reproduced by keeping only the lowest open string modes. These 
phase shifts corresponds to an eikonal approximation where the gravitons move along 
infinite straight lines at large impact parameter. Similar calculations have later successfully 
extended this to processes involving higher branes (as in |5|,||[7|]) and spin effects ||, as well 
as processes involving longitudinal momentum transfer ||. The Matrix theory calculation 
of the graviton scattering process has lately been pushed to two loops fK3,|llJ , giving further 



evidence for the finite N conjecture. Recently Seiberg |E| has given arguments for why 



the Matrix theory conjecture is correct, and this was further examined in the context of 



graviton scattering in |13 



From the infinite time phase shift calculation in Matrix theory one derives an effective 
potential between gravitons of relative velocity v and relative distance r. This potential 
is a double expansion in i> 2 /r 4 and 1/r 3 (the latter is the loop expansion), starting out as 
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ff = I6 K V + ° ' " ' V + 4096* ' V + ° <K V 
+ 0-^ 3 ^ + § K ^ + O(- t »^) (1-1) 

+ 0(" 14/3 ^). 



where we have indicated how each term scales with the 11-dimensional gravitational cou- 



pling, otherwise (when setting k = 1 above) the units are as in |10[. So far, the two first 
terms along the diagonal of this double expansion, with integer powers of k, have been 
found to coincide with supergravity [II |. It has been argued that all the terms along the 



diagonal will agree with classical supergravity and that all terms to the left of the diagonal 
should vanish |Il,|IRO], while other non-zero terms in this expansion should come from 



higher-derivative terms of the supergravity effective actions |13 . 

These infinite time straight line phase shift calculations remain one of very few quan- 
tities one can actually calculate directly in Matrix theory. It is important to try to extend 
this repertoire as much as possible. In this paper we will consider a finite time version of 
this calculation, and from the resulting "phase shift" we can read off an effective potential 

locally, rather than integrated along an infinite path. We will only work to one loop, and 

i <i 4 

there we find again the y| 77 potential of supergravity to leading order, for arbitrary time 
intervals. This provides a stronger equivalence between Matrix theory and supergravity 
than is demonstrated from the infinite-time the finite time amplitude contains 

more information. Of course, the main conjecture of Jl| concerns the 5-matrix, so this re- 
sult suggests that their conjecture could perhaps be strengthened. However, higher order 
terms in the potential will be modified compared to the infinite-time case, by interesting 
terms that integrate to zero along the infinite line, thus showing that the complete match- 
ing of the leading term is non-trivial. Working on a finite time interval means that we have 
to be careful about what boundary conditions are put on the high energy modes that are 
integrated out. We will consider the most natural boundary condition for comparing to 
supergravity, and also investigate how the answer varies with other choices. We find that 
the leading term is robust towards changing the boundary conditions, while the subleading 
terms are more sensitive. In the large iV limit proposed in , these terms are also sublead- 
ing in powers of N, but they may be interesting in Susskind's finite N conjecture ||. Note 
that if the Planck scale lp is the only relevant scale in supergravity, then since k^{ 3 ~ lp 
these terms seem to indicate corrections to the supergravity action dimensionally going 
like R + R 5 / 2 + ... rather than R + R 2 + .... 

This paper is organized as follows: In sect. 2, we set up the basic problem. Sect. 3 
contains our calculations, and sect. 4 has some concluding remarks. 



2. Comparing Matrix theory and supergravity 

Graviton scattering in Matrix theory is described by U(N\ + N2) supersymmetric 
quantum mechanics, where N\/R and N2/R are the longitudinal lightlike momenta of the 
gravitons, R being the radius of the compact null-direction. We will take Ni = N2 = 1 
here, as usual the leading order factors of Ni can be easily reinstated afterwards. The 
bosonic part of the U(2) action is given by ]I7J 

S = TrJ dt(±(D t X i ) 2 + ~[X i ,Xi}A, (2.1) 

where for simplicity we suppress dependences on R and the 11-dimensional Planck scale 



see 



0). We defined D t X % = d t X i + [A,X% with A, X % being U{2) matrices which we 



can decompose as 



X* = % - (X* I + I> a ) . (2.2) 



The o" a , a = 1 ... 3, are the Pauli matrices. There are also corresponding fermionic fields 
tjj. The Xq fields describe the center of mass motion, and they (together with Aq and i/jq) 
decouple from the rest and will be ignored from now on. We interpret X$ as the relative 
coordinate of the gravitons, while the corresponding gauge field As can in principle be 
gauged away (in background field gauge, which we employ here, it is set to a constant). 
The fermionic ^3 fields incorporate the spin and polarization of the (super-)gravitons. 
The off-diagonal 1,2 modes do not have such an intuitive description in the long distance 
supergravity, in fact they can be thought of as arising from open string considerations valid 
only at short distances, and they then represent the unexcited modes of these open strings. 
When the gravitons are far part, the off-diagonal modes are all very heavy, with masses 
roughly proportional to the distance. One can then imagine integrating these modes out 
in a Born-Oppenheimer type of approach, as (non-local) internal degrees of freedom. This 
leaves a theory only involving the diagonal modes that are interpreted as the positions of 
supergravitons. This separation of light and heavy modes is only exact in the limit of large 
<X\> though, and a more systematic understanding of this is needed, especially at higher 
loop orders. 

When integrating out the off-diagonal modes, we need to supply boundary conditions 
on the fields. In previous infinite time phase shift calculations, the fields are taken to 
vanish at ±00, but for finite time intervals we should be more careful. The off-diagonal 
modes are at one loop level described by harmonic oscillators, with frequency proportional 
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to their mass, and specifying boundary conditions can be done by specifying the initial and 
final states of the harmonic oscillators. These harmonic oscillator modes do not appear 
at all in supergravity, so exciting them out of their ground state would seem to take us 
out of the realm of supergravity. For comparing finite time amplitudes, the most natural 
assumption therefore appears to be to put these modes in their ground state for both the 
initial and final states. We will also consider other choices of boundary conditions. 

If we also specify boundary conditions for X 3 and integrate it out, we will find the 
amplitude for a finite-time graviton propagation. We can write this schematically as (we 
now work in Euclidean time) 

(Xi s ,T s \X l 3J ,T f ) = Af VX a e~ s{Xa) =M VX 3 e~ s ^ X3) = e~ So+s . 

J X 3 (T s )=X3 jS J X 3iS 

(2.3) 

Here Sq is the action of the classical trajectory implementing the boundary conditions, 
and in an abuse of language we call S the finite time phase shift and try to relate it to 
an effective potential through 5 = — f V e R. Since we are not considering any sort of spin 
effects, ips (as well as A 3 , the ghost C 3 and X 3 once the boundary conditions have been 
implemented) just go along for the ride as massless free fields here at one-loop level. 

In supergravity the most convenient way to do a similar computation, is to give one 
of the gravitons a large longitudinal momentum such that it can be treated as a classical 
source for the other graviton. One can then derive an action for the probe graviton moving 
in the field of this source, 5'sg(^3)- There are then two ways to view the comparison 
between Matrix theory and supergravity. On one hand, we can directly compare the 
effective actions of the light modes S e s(Xl) and Ssg(^3) f° r various paths X^(t) and see 
how they match — our computations in this paper are equivalent to doing this for a finite 
line segment. This seems to be the cleanest way of interpreting our results. Alternatively 
we can evaluate the finite time amplitude in supergravity, which, at the semiclassical level, 
is given by integrating the action along the geodesic connecting the two space time points. 
At large separation of the gravitons, the leading order amplitude can be found using an 
eikonal approximation, replacing the geodesic by a constant velocity straight line. 

We will choose coordinates such that the straight line connecting the two space time 
points is given by 

Xl (t) = vt, 

X!, (t) = b, (2.4) 
X^(t) = 0. 
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The supergravity action for the probe graviton is given in and we find that 
the supergravity action evaluated for a finite time line segment (alternatively, an eikonal 
approximation to the finite time phase shift) is given by 



15 v 4 



+ 0(- 



.14 



(2.5) 



16r 7 (t) 

where r(i) = vW+F. The 0( p-j) term we don't expect to see until two loops in Matrix 



theory and will not concern us here. Note that, as argued in ||18|| , the corrections to the 
eikonal approximation are of order v 2 /r 7 and thus fall along the diagonal in ( |1 . 1| ) and will 
not be of concern to us here, either. 



3. Finite-time calculation in Matrix theory 



On the Matrix theory side, to evaluate (|2.3|) we implement the boundary conditions 
on XI by a suitable solution to the equations of motion. These solutions are straight lines, 
so we expand about X\ Q (t) in (|2.4| ). The path integral can then be evaluated as in ||, at 
one loop it yields a product of determinants. The phase shift is given by 



5 = ln[det" 6 (-a t 2 +r 2 ) det _1 (-d t 2 +r 2 +2i;) detr^-d^H 
x det 4 (-<9 2 +r 2 +t0 det 4 (-a 2 +r 2 -t;)]. 



-2v) 



(3.1) 



Now these determinants represents the evolution of a time dependent harmonic superoscil- 
lator from a state \ijj s ) at time T s to a state \i^f) at time Tf. The boundary conditions 
we put on the fermionic part should be determined by supersymmetry from the boundary 
conditions on the bosonic part. We here choose to represent each determinant by a sepa- 
rate bosonic harmonic oscillator of the appropriate frequency, evolving from corresponding 
initial states to final states. This prescription preserves a symmetry between all the modes, 
and we believe it gives equivalent results to treating the full superoscillator. We write this 
as 



det~^{-df + u 2 (t)) = (^/|Texp 



H(t)dt 



M / dipfdipsipUipf^sip. 



T>cp(t) exp 



L(t)dt 



N I dffd<f s i/jf(f f )'il; s ((f s )exp 



deV(-d 2 + c 2 (t)), 



(3.2) 
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where 



H 
L 



(p 2 + u 2 ip 2 ) , p 

(v? 2 + ^V 2 ), 



-id. 



(3.3) 



[-d 2 + u 2 {t)\ <p (t) = 0, (f (T a ) = ip 8 , (po(T f ) = ip } . 

By deto we denote the determinant evaluated on the space of functions that vanish at both 
endpoints. 

In order to extract the leading order supergravity potential, we actually need to do 
very little work, the simplest adiabatic approximation is sufficient. In this approximation, 
we assume the harmonic oscillator to stay in its ground state throughout, and we find the 
estimate 



det"2(-d 2 + u 2 (t)) « exp 



> 1 
-u{t)dt 



(3.4) 



Using uj 2 {t) = v 2 t 2 + b 2 + av, inserting this into ( |3.1| ), and expanding in powers of v, we 
get the phase shift 



>ad 



(it 



15 v 4 315 v 6 27027 v 8 

16 + 128 + 4096 



+ ... 



(3.5) 



We see that the leading term at large distance agrees with the supergravity result (|2.5|), 
while the higher terms do not even agree with the infinite time phase shift calculation. 
This way of deriving the potential energy between the gravitons is tantamount to summing 
up the zero point energy of the harmonic oscillators representing the off-diagonal modes, 
which is how various potential energies were computed for instance in 0. We see that at 
subleading order there is a distinction between these two methods of extracting an effective 
potential. 

We will now embark on a more systematic evaluation of the finite time amplitude. 
We will use two approaches, valid, roughly speaking, at long and short time intervals 
respectively. 

For long time intervals, we use the last line of (|3.2|) . The classical trajectory ^o(^) is 



a parabolic cylinder function ||19[| . We can systematically solve for it using a WKB type 
expansion (equivalent to the Darwin expansion in jn|). To this end, write 



<p (t) = exp 



/ Y^n n - X f n {t')dt' 

JT, n 

a n=U 



(3.6) 
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and consider u(t) to be of order 1/ft. Solving the equation for (p in (|3.3|) order by order 
in h (and then setting ft = 1), we find the first few /'s 



fo = u 
fi = - 



vH 



2a; 2 ' 

f 2 = ^(2u; 2 -5v 2 t 2 ), ( 3. 7) 
/3 = ^(3^-5^), 



4 



/ 4 = — (—76a; + 884wVr - 1105vV). 
128w 

There is also another solution obtained by f^n — > —f2n for all n. We see that 

/ v \ n 

fn~u J g n (vt/u), (3.8) 

where g n only contains non-negative powers of its argument. Thus this is an expansion in 
v/uj 2 which is what we want. If we define 

Ait) = J2f2n(t), Bit) =J2hn+i(t), (3.9) 
the general solution can be written 

Mt) = cJ> +A) + cJ>~ A) . (3.10) 
Solving the boundary conditions, we can evaluate 



ip (T s ) = Lp s B s - cp s A s coth j A + iPf J^£r A 



f A I B 

<p (Tf) = <PfB f + <p f A f coth J A- °^r A ■ 



(3-11) 



Here A s = A(T S ) and so on. It is now straightforward to evaluate the integral over f s j in 
(|3.2|) . The ground state wave function is 

/ U) \ V 4 1 2 

iKv) = (-) e-* u * , (3.12) 

and the integral is gaussian. The remaining determinant deto(— d 2 +w 2 (t)) can be evaluated 
using the method explained in |16| , p. 340. The prescription is to solve the equation 

(-d 2 +u 2 (t)) X (t) = ] X (Ts) = 0,x(Ts) = l, (3.13) 
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then 

det Q (-d 2 +u 2 (t))=M' X (T f ). (3.14) 
We already solved this equation above. Solving for the boundary conditions, we find 

eJ 



x(T f ) 



A* 



sinh / A 



1 



^A~A 



sinh / A 



f 



(3.15) 



using B = —A/ (2 A) which is easily derived from the equation. 

We now restrict to the regime where J 10 ^> 1, we can then replace all hyperbolic 
functions by exponentials up to errors of order exp(— 2 f u). Doing the gaussian integral, 
we find 



det(-df + u 2 (t)) =M'' 



1 + 



2(x! a 



1 + 



A'+B f 



2lo 



f 



A' 



A' ' 
1 + - L 



exp 



J(u + A') [l + 0(e- 2 I l 



M" exp 



f , 1 
ujdt + - 



> vH 2 



-dt 



(3.16) 

Here we defined A' = A — fo and B' = B — f\. The result takes a simple symmetric form, 
every second term in the expansion can be written as an integral, while the others have 
equal contributions from each endpoint (the first of these actually cancels completely). 
This latter type of term could in principle be removed by a suitable normalization of the 
initial and final wave functions we use, however as far as we can see, there is no rationale 
for doing that. There are no ambiguous time-dependent normalization factors floating 
around here, as we see from comparison with the short-time calculation below. When we 
write 0(v 3 /w 6 ), it is meant up to positive powers of vt/u. We can now plug this into (|3.1J ) 
and expand in powers of v, this gives the phase shift 



T ' 15 , 
T 16 r' 

J s 

T f 429 v 8 



f 315 v 



6 



45 



T 2 



128 r 



-(l-ll(f.^+-^ + 



(3.17) 



-I / ^\63 + 1500(r-v) 2 - 12070(r-v) 4 ]dt + O(v 8 /r 16 ) + O(e~ 2 f r ) 

4096 r Lt3 



Here f = f/r. This result reduces to the first line of ( |1 . 1| ) in the limit T s — > — oo, Tf — >• oo 
as it should, but we see the interesting feature that it differs for finite time intervals. The 
leading order potential still matches supergravity perfectly. The subleading potential is 
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different from both the instanteneous zero-point energy potential and the potential derived 
from the infinite time phase shift — in a sense it interpolates between these two. There 
are also terms depending on the end points additively and thus cannot be written as an 
integrated potential. The interpretation of these terms is not clear. Note that the first 
subleading term is a simple total derivative, 

^{l-ll{r-v?) = j t ^-v (3.18) 

for the straight-line trajectory. It seems unlikely that the corresponding term for a curved 
trajectory would be a total derivative [^D| . 

It is interesting to see how the result changes for different choices of boundary con- 
ditions on the off-diagonal fields. If we choose boundary conditions such that the fields 
vanish at the end points, the determinants are just given by deto, and we find the phase 
shift 

— — (l-ll(f.0) 2 )^ + .... (3.19) 




This also reduces to (|1.1| ) in the infinite time limit, as expected, and the leading potential is 
unchanged, but the subleading terms have changed. We can also imagine picking boundary 
conditions such that we have the ground state initially while summing over all final states, 
this corresponds to putting ifjf(Lpf) = 1 in (|3.2|). Also in this case we find the same leading 
order potential, while the higher terms are now asymmetric in Tf and T s and cannot be 
written in a particularly illuminating form. 

The approach above breaks down for AT = Tf — T s very small. We want to check 
further what the Matrix theory predicts in this regime. In the infinite-momentum frame, it 
is not clear to us whether physics at short time scales as described by Matrix theory should 
agree with physics at short time scales as described by supergravity. Since supergravity is 
an effective low-energy description of M-theory, one might expect Matrix theory to behave 
differently from supergravity in this regime. To deal with short time-intervals, we use the 
first representation of the determinant in (|3.2|). We write 



H = H s + H p , 

H s = ^(p 2 +w s V), (3.20) 

n v = \v\?-T^\ 



and treat H p as a time-dependent perturbation. This expansion should be useful for 
v 2 (T 2 — T 2 ) small, in particular u s AT can be allowed large thus giving an overlap in 
the region of validity with the expansion considered above. We expand the time ordered 
product as 

det-i (-d 2 + u\t)) = (0/| e"^ AT |0.) - / ' (0/| e'^^ H p {t)e- H ^- T ^ %) + .... 

(3.21) 

Here \n s ) denotes the energy eigenstates of the harmonic oscillator at time T s . To evaluate 
this expansion, it is convenient to expand |0/) in terms of the \n s ). Representing ip 2 in 
terms of creation and annihilation operators it is then a matter of straightforward, but 
tedious algebra to evaluate the determinant. The result for short time intervals is 



det(-d 2 +w 2 (t)) = exp 



Ujdt + lA^ AT 2 + l 7 - 



8^ 



24ivt 



v 2 T 2 , 

-2<X t 3 - ) AT 3 + C(AT 4 ) 

(3.22) 



We have checked the consistency between this expansion and the WKB expansion above 
in the region where v 2 T s AT < 1, uAT > 1 (we checked all terms involving t> 4 ). From 
the determinant we calculate the short time expansion of the phase shift, 



AT 



15 v 4 315 v 6 

16 rl + 128 ^l 1 + 



+ 



-15- 



,6 



,10 



105 ^ + 

s 



+ AT 2 



105 v l 
+ C(AT 3 ) 



3465 v 1 
256 



+ 



[f s 



(3.23) 



We see that even at short time we reproduce the supergravity result to leading order, 
the piece linear in AT is nothing but the simple adiabatic result from above, while we 
see deviations from this approximation at order AT 2 . These deviations cannot be easily 
written in form of a potential, which would have required 5 to take the form 



5 = 



T a +AT 



V eS (t)dt = -V eS (T s )AT - ^ eff (T s )AT 2 + . . . 



(3.24) 



The new terms appearing at order AT 2 are indicative of terms in the effective Lagrangian 
that are not expressible as a single local time integral. 

If we consider the short time behavior for the case of vanishing boundary conditions 
on the off-diagonal fields, we find a very different result. We can evaluate the short-time 
determinant by expanding %(t) in a Taylor series, and solve (|3.13| ) order by order in AT. 
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The first terms in the determinant will then be proportional to cu 2 , cu 4 and so on, and these 
terms all cancel when multiplying the determinants for the phase shift. In fact the phase 
shift vanishes all the way up to AT 8 and we find 

S vv = -^-AT 8 + 0(AT 10 ). (3.25) 

This strange result probably shows that this is a particularly unwise choice of boundary 
conditions at short time intervals. 



4. Discussion 

When considering finite time amplitudes in Matrix theory, one must be careful in 
specifying boundary conditions on the off-diagonal modes that have no counterpart in 
supergravity. The most natural choice seems to be to demand these modes to be in their 
ground state. In this case we derive an effective potential between the gravitons that 
differs from the one obtained from the infinite time calculations. The leading order term, 
which is the one reproduced by supergravity, still matches perfectly, providing a functional 
agreement between supergravity and Matrix theory. In comparing these Matrix theory and 
supergravity calculations, one expects the low energy considerations to be valid only away 
from very short time intervals. We find that for long, but finite time intervals, there is a 
correction to the effective potential, depending on the angle between r and v, compared 
to the one read off from the infinite time phase shifts. There are also new terms that 
are not expressible in terms of a potential, whose role is not so clear. In particular the 
first subleading term in the potential, of order v 6 /r 11 , is non-zero. Note that this term 
dominates the w 6 /r 14 -term (which has been matched between two loop Matrix theory and 
supergravity) at long distances. Since supergravity should be corrected by M-theory effects, 
it would be interesting to interpret the term we found as such an M-theory correction. 

This type of calculation should be straightforward to extend to include spin effects and 
other objects than gravitons. We expect the same features to show up also then, namely 
that the leading term matches supergravity as in the infinite time case, with corrections 
to the subleading terms. It would also be interesting to consider extensions to higher loop 
calculations. 
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